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SUMMARY 


An application of a simple numerical technique which allows for the 
rapid construction of orthogonal coordinate systems about axisymmetric blunt 
bodies is presented. This technique can generate orthogonal meshes which have 
unequally spaced points in two directions. Relations are given for the numer- 
ical generation of the metric coefficients. 

Body shapes ranging from simple analytical bodies to complex reverse cur- 
vature bodies are presented together with their orthogonal coordinate systems. 
The relatively good accuracy of the technique is shown in tabular data describ- 
ing coordinate line slopes and metric coefficients. The "predictor-corrector" 
numerical method used to generate these results is both simple in concept and 
easy to program, so that the application of the technique should be broader 
than the results presented. 


INTRODUCTION 

The introduction of high-speed computers with extensive memory capabilities 
has greatly enhanced the feasibility of solving the full Navier-Stokes equations 
for fluid flow over complex geometric shapes. One of the primary problems hin- 
dering the development of computational fluid dynamics for complex geometries 
has been the difficulty of generating the finite-difference mesh. A great deal 
of effort has been expended to develop coordinate transformations and/or mesh 
generators for varying degrees of geometric complexity. (See refs. 1 to 14.) 
These techniques range from the very simple, as in reference 1 , to the mathemat- 
ically elegant, as in reference 14. However, for flow over blunt entry bodies, 
it is imperative that the boundary conditions on the body be represented as 
accurately as possible by the finite-difference approximations to the Navier- 
Stokes equations, because the regions adjacent to solid surfaces generally domi- 
nate the character of the remaining flow field. The viscous shear stress forces 
on solid bodies depend directly on the large gradients (in the direction normal 
to the solid surface) . Accurate skin-friction coefficients require that these 
large gradients be resolved accurately. Therefore, almost all numerical solu- 
tions to the Navier-Stokes equations to date have used body geometries condu- 
cive to use with natural or nearly orthogonal coordinate systems. 

In the natural coordinate system, the body surface itself forms one bound- 
ary. That is, the body contour coincides with a constant coordinate line. Some 
simple examples of this approach are a cylindrical coordinate system to describe 
flow over a cylinder, a spherical coordinate system to describe flow over a 

sphere, and a parabolic coordinate system to describe the flow over a paraboloi- 

dal body. In each of these coordinate systems, the normal coordinates intersect 
the body orthogonally, thus simplifying the boundary conditions. There is lit- 
tle difficulty in compressing the mesh near the body because the computational 
mesh system is composed of lines parallel to the body which can be concentrated 

as close to the body as desired. There is, however, one rather severe restric- 



tion on the natural coordinate system; that is, the body must have an analytic 
shape. Unfortunately, most realistic entry bodies bear little resemblance to 
the limited number of natural coordinate systems available. 

Another option is similar to the natural coordinate system in that the 
body surface becomes one coordinate line in the system. This is called the 
body-oriented coordinate system. In this system the coordinates of a point are 
determined by the distance along a body surface measured from the axis of sym- 
metry and the distance along a normal to the body. This type of system has 
been used to describe the flow over the forebody portions of blunt entry bodies. 
(See refs. 15 to 17 for examples.) The major problem with this system is that 
any discontinuity in slope on the body makes it impossible to describe the com- 
plete flow field. The body-oriented system, therefore, cannot handle the body 
shape with a concavity produced by heat-shield ablation in a severe heating 
environment. 

Conformal mapping gives a well-developed and accurate alternate method 
of generating orthogonal coordinate systems for complicated geometric shapes. 
Reference 1 4 provides a recent look at some of the latest and most sophisticated 
techniques. Conformal mapping methods, however, are mathematically complicated 
and require multiple transformation steps leading to a loss of physical reality 
in the computational plane. Such complications make finite-difference mesh 
setups difficult. 

Reference 12 presents a nearly orthogonal coordinate generator, somewhat 
similar to conformal mapping in concept, which has been used in two-dimensional 
and axisymmetric calculations; however, the technique is not easily applied 
and requires a great deal of user familiarity to execute successfully. One 
relatively easy analytic technique (ref. 1 0) has produced excellent results 
(ref. 18), but the technique is limited by the necessity of specifying the 
body shape by an analytical trigonometric series. 

This paper presents a more general orthogonal coordinate system which can 
represent a wide range of body geometries. This technique is based on the 
rather simple but accurate numerical method developed in reference 19. The 
body can be represented by a series of discrete (but continuous) points rather 
than by an analytical approximation. In the transformed computational plane, 
the region of interest is rectangular, with the body surface being a coordinate 
line. This representation combines the advantages of the natural coordinate 
system and the generalized orthogonal coordinate system. In addition, generat- 
ing the coordinates in the physical plane simplifies the finite-difference mesh 
setup. 


SYMBOLS 

hi »h 2 ,h 3 metric coefficients for transformed coordinates 

hi,j metric coefficient 

i tangential direction index 


2 



J 

j 

Kc 

N 

r Z 

r o 

r s 

U 1 ,U 2 ,U 2 
X 1 ,X 2 ,X 2 

x,y,z 




0 

e o 

e i 

p 


total number of mesh points in normal direction 

normal direction index 

unequal spacing parameter 

unequal spacing coordinate parameter 

local distance between body surface and outer boundary 

radius of outer circular boundary 

radius of body surface from its geometric center 

generalized transformed coordinates 

generalized coordinates 

Cartesian coordinates 

transformed orthogonal coordinates 

internal body angle 

outer circular boundary angle 

local angle for r^ 

local radius from axis of symmetry 


Subscripts: 

o outer boundary related quantities 

s body surface quantities 


DESCRIPTION OF METHOD 


For the mesh generation about an axisymmetric body, the origin of the 
x,y coordinate system is taken as the approximate center of the body. Once the 
body is described in terms of x,y points, then the surface distance E, can 
be calculated in a generalized manner by defining £ as 0 on the -y ordinate 
(see fig. 1 ) and 0 as increasing positively in the E, direction. Then E, 
is given by 


5 ( 0 ) 




d0 
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where 


r s = (x s 2 + y s 2 ) 1/2 
6 = sin -1 (x s /r s ) 


Once K is defined, the next step is to construct ri level lines between 
the body and a suitable outer boundary. As pointed out in reference 11, a 
circular outer boundary simplifies the specification of the outer boundary 
conditions for fluid dynamic calculations. Most of the examples presented 
in this paper use a circular outer boundary. However, the outer boundary 
may be specified as any reasonable shape; and for those cases where only a 
partial field is to be described, one can use a discrete shock wave as the 
outer boundary. The circular boundary can be easily specified by a polar 
coordinate system where 


x = r Q cos 0 O 


y = r Q sin 0 O 


In this system, r Q is the radius of the outer boundary and 0 O is usually 
taken in equal increments around the outer boundary. On the outer boundary, 
h = 1 while on the body ri = 0. The level lines between the boundary and the 
body can be constructed along straight lines connecting corresponding points 
on the body and circular boundary. Note that the mesh points on the outer cir- 
cular boundary are not the final mesh points, but initial values used only to 
set up the level lines. The spacing of the level lines can be unequal and the 
unequal spacing relation of reference 20 can be easily applied: 


^Nj/AN _ 1 

n i = : 

K^-l 

where Nj = (j - 1)AN and Nj = 1 , with AN = 1/(J - 1) and K c being the 
spacing parameter (generally 1 < K c < 2) . The larger the spacing parameter 
K c , the more unequal the spacing. Frcxn the above relation, the level lines 
between all corresponding points on the body and circular boundary can be 
calculated. The relationships for the corresponding values of x,y can be 
obtained (see fig. 2 geometrical schematic) from 


2 2 

r 2 , i = (*n=i - x n=o> i + (Yn=i - Yn=o). 


1/2 
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where 


x i,j = x i,n=0 + ^j r Z,i) cos < e l) 
Yi,j = yi f n=0 + sin < e l) 


01 = sin" 1 [(y n=1 - y n =0)Aj,i] 

Figure 3 shows the level lines constructed in this manner for an ellipse with 
K c = 1.01, which gives nearly equal spacing. 

Once the level lines have been determined, the normal lines are con- 
structed so that an orthogonal system is defined. The approach to the con- 
struction of the normal lines is the one given in reference 18 which uses a 
simple "predictor-corrector" technique analogous to the trapezoidal integration 
technique of numerical integration. In this technique, the solution is first 
predicted from the level line at a known point by using the Euler method. Once 
the predicted point on the next level line is obtained, then the slope at that 
point is calculated and a new predicted point is obtained from this slope. The 
actual solution is then a combination of these two solutions; i.e., the final 
x and y values are an average of the predicted and corrected values. This 
procedure is illustrated in figure 4. The solution then proceeds point by 
point along a level line until all normals on that level have been constructed. 
Then the solution proceeds to the next level and the process is continued until 
the outer boundary is reached. 


EXAMPLES OF COMPUTED COORDINATE SYSTEMS 

Figures 5 to 7 show examples of the orthogonal coordinate systems con- 
structed about three ellipses with increasing degrees of ellipticity. These 
solutions are for equal spacing both in the £ and U directions. Notice 
the divergence of the coordinate lines in the regions of large curvature. 

This effect of curvature can be easily remedied through unequal spacing in 
the £ direction as shown in figures- 8 and 9. If unequal spacing is used 
in the T) direction also, then a smoothly varying coordinate system with no 
region of rapid divergence of the coordinate lines results. (See fig. 10.) 

A somewhat more realistic body is presented in figure 11. This body rep- 
resents a typical planetary entry shape, and the coordinate system produced is 
quite good even though very little unequal spacing is used in the £ direc- 
tion. Figure 12 shows the same body, but with unequal (K c = 1.08) spacing in 
the ri direction. Again, as in figure 10, the results are quite good. 

To demonstrate the use of the coordinate generator on bodies with reverse 
curvature, the body shown in figure 13 was arbitrarily chosen. The reverse 
curvature poses no special pr oblem; however, the normal coordinate lines tend 
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to converge. The effect is highly unequal spacing in the regions far from the 
body. The amount of reverse curvature here is relatively small; nevertheless, 
problems will probably occur on deeper cavity shapes. Figure 14 shows the same 
body with unequal (K c = 1.08) spacing in the ri direction. 

A natural compression of n lines occurs when the geometric center of the 
outer boundary is displaced in relation to the geometric center of the body, as 
shown in figures 15 and 16. For supersonic blunt body flows, the shifting of 
the outer boundary is highly desirable, since the shock wave is close to the 
body in the forward region and a compression of coordinate lines is necessary. 
The stretching of the level lines in the base is not detrimental, since the nor- 
mal gradients in this region are not as strong as on the forebody. 

All the above cases were run with 31 n level lines and 100 £ normal 
lines giving 3000 points for the present analysis. Run times on the CYBER 175 
computer averaged approximately 5 seconds, including compilation. On more com- 
plicated geometries, more points may be necessary. For cases where there is a 
symmetry plane in the field, only half as many points need to be used. 

In addition to the whole body solutions given above, the coordinate genera- 
tor can be used on partial bodies such as the forebody region on a blunt body. 
Figures 17 and 18 give solutions for partial bodies, with a representation of an 
ablated entry body with reverse curvature in the stagnation region shown in fig- 
ure 18. This technique is sufficiently fast to make its use with an interactive 
graphics terminal both feasible and desirable. 


COORDINATE METRICS 

Since the intended use of the orthogonal coordinates generated using the 
present technique is for axisymmetric bodies, the following coordinate system 
(see fig. 19) will be used to generate the metrics, where the nomenclature fol- 
lows that of reference 21 : 


c 
— 1 

II 

Z3 

x 1 = p 

COS 

4> 

u 2 = K 

x 2 = p 

sin 

<t> 

u 2 = 4> 

X 2 = z 



where p = P (5,n) ? 

N 

II 

N 

?M 

and 

the metric coefficients are given by 


if 3 


3x^ 3x^ 3x^ 0 X 2 0 x 3 3 x 3 

■ — “ + 1 1 ■ + — — 

3u* 3u7 3u* 3u7 3u* 3u7 
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For an orthogonal system, the metric coefficients 


h l , 2» h 2,lf h l,3* h 3,l t h 2,3» h 3,2 


all have to be zero, leaving only the three familiar coefficients h-j ,h 2 ,h 3 , 
which are 


h] = hl,l h2 = h 2/ 2 h 3 = h 3,3 


When the derivatives in the metric coefficient relation are taken, the following 
metric coefficient expressions are obtained: 



h 1 ,2 
h l ,3 
h 2, 2 


3p 3p 3z 3z 

3rj 3^ + in i? 


h 3,l ° 0 



h 2, 3 = h 3,2 = 0 
h 3 , 3 = P 2 


With the choice of the present coordinate system, two of the three necessary 
orthogonality conditions are identically satisfied, leaving only the h-j 2 
coefficient to be evaluated numerically to determine the accuracy of the f orthog- 
onal coordinate generator. Table I presents results for the metric coefficients 
evaluated with second-order accurate derivative relations. Generally, h-j 2 is 
relatively small compared with the other coefficients in spite of being a pro- 
duct of derivatives which amplifies errors. One has reasonable justification 
for ignoring the hi 2 coefficient and assuming orthogonality; thus only the 
three coefficients , ^2,2t an ^ ^3,3 are left. Orthogonality is also 
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indicated by comparison of the slopes of the level and normal lines at their 
intersection. Table II gives some typical results which indicate that the 
normal-line construction procedure does, in fact, give reasonably good results. 

To test the effect of mesh size on the metric coefficients, a 1.5:1 ellipse 
was chosen with equal spacing along and normal to the body. The hi i metric 
coefficient on the body surface is given in table III for a range of mesh sizes. 
Because of symmetry, only the results for the first quadrant are presented. As 
expected, the results show that the metric coefficient does vary with mesh size; 
however, for the range of mesh points chosen, the variation is relatively small. 

The only significant drawback to the numerical coordinate generator scheme 
is that the metric coefficients also must be numerically generated and thus 
stored at every point for future use. Since the metric values must be stored 
at each point, some computer storage limitations could result for very large or 
complex geometries. 


CONCLUDING REMARKS 

The results of the application of a simple "predictor-corrector" numeri- 
cal technique which allowed rapid construction of orthogonal coordinate sys- 
tems about axisymmetric blunt bodies have been shown. This technique generated 
orthogonal meshes which have unequally spaced points in two directions. The 
technique is general enough to allow for coordinate generation about bodies 
ranging from simple analytic shapes to more complex bodies with reverse cur- 
vature. The relatively good accuracy of the technique was demonstrated through 
the tabular data presented on coordinate line slopes and metric coefficients. 

The predictor-corrector numerical technique used to generate results given 
in this paper is both simple in concept and easy to program, so that the use of 
the method should be broader than the limited application presented here. In 
fact, the method should lend itself to use with interactive graphic systems. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 17, 1979 
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TABLE I.- METRIC COEFFICIENTS 


(a) Ellipse; £ = 1.34 


n 

h l,l 

hi, 2 

h 2,2 

h 3,3 

o. 

.29743E+02 

-.52567E-02 

.99942E+00 

,792888+00 

.28748E-01 

,313648+02 

-.377428-02 

, 10479E + 01 

.10957E+01 

.577848-01 

. 325 1 6E+02 

-.453428-02 

, 1 108SE + 01 

,145738+01 

,871108-01 

.33369E+02 

-.310528-02 

,118108+01 

.18789E+01 

,116738+00 

.34036E+02 

-.191228-02 

.12637E+01 

.23620E+01 

,166648+00 

, 34531E+0? 

-.772458-03 

.13567E+01 

.29078E+01 

,176868+00 

,349218+02 

-.128448-02 

.14601E+01 

.351758+01 

,207388+00 

.352358+02 

-.45S89E-03 

, 15733E + 0 1 

,419238+01 

.238208+00 

.35468E+02 

, 28407E-02 

,169628+01 

.493348+01 

,269338+00 

.35648E+02 

.235478-02 

.1B298E+01 

.57418E+01 

,300778+00 

.35785E+02 

.288118-02 

, 19744E + 01 

.66185E+01 

.332528+00 

,358878+02 

,301918-02 

.21301E+01 

.75647E+01 

,364608+00 

,359618+02 

.313528-02 

,229788+01 

.85815E+01 

,396998+00 

.36013E+02 

.326428-02 

, 24778E + 0 1 

.96698E+01 

,429718+00 

,360«5E+02 

.33919E-02 

, 267 07E + 0 1 

.108318+02 

.462768+00 

,36 »62E+02 

,351828-02 

.287718+01 

.12066E+02 

,496138+00 

, 36066E+02 

,365048-02 

.30976E+01 

, 1 3375E+02 

,529848+00 

,360588+02 

, 37923E-02 

.33329E+01 

,147618+02 

,563898+00 

.36041E+02 

.394608-02 

.35834E+01 

.16224E+02 

.59827E+00 

,360 1 5E+02 

'.41 149E-02 

, 38500E + 0 1 

.17765E+02 

,633008+00 

,359838+02 

.43034E-02 

, 4 1 333E + 0 1 

.19385E+02 

.668088+00 

, 35944E+02 

.450778-02 

, 44339E + 0 1 

.21085E+02 

,703518+00 

.35900E+02 

,472178-02 

, 47526E+0 1 

.22867E+02 

,739308+00 

, 3585 1 E + 02 

.552788-02 

,508998+01 

.24731E+02 

,775448+00 

.35798E+02 

,574408-02 

,544638+01 

.26679E+02 

.811948+00 

,35741 E+02 

,5951 18-02 

.S8229E+01 

.287118+02 

,848818+00 

.356818+02 

,614508-02 

.62204E+01 

.30829E+02 

.88604E+00 

, 3561 8E+02 

.632108-02 

, 66397E + 0 1 

.33034E+02 

.923658+00 

, 35553E+02 

, 6 1 864E-02 

.70B16E+01 

.35326E+02 

,961648+00 

, 35485E+02 

-.243488-02 

.75479E+01 

.377088+02 

.100008+01 

.35415E+02 

.43529E-02 

.80379E+01 

.40179E+02 


11 










TABLE I.- Continued 


(b) Ellipse; K = 3.72 


n 

h l,l 

CN 

-C 

h 2,2 

h 3, 3 

^ • 

,137396+02 

-.232596-02 

,100156+01 

.1 1855E + 00 

.287U8E-01 

,196036+02 

-.989176-02 

, 1 2339E+0 1 

,180996+00 

,577896-01 

,157836+02 

-.150506-01 

, 1955 1 E + 0 1 

.25929E+00 

,871106-01 

, 16989E+02 

-.175796-01 

,167326+01 

,356996+00 

,116736+00 

.18J03E+02 

-.171356-01 

.18931E+01 

.97269E+00 

, 146646+00 

,191666+02 

-.187956-01 

.211786+01 

,608916+00 

, 1768bE+00 

.20099E+02 

-.189996-01 

,23988E+01 

,763966+00 

,207386+00 

,208906+02 

-.219696-01 

, 2589 1 E + 0 1 

.93897E+00 

, 238206 + 00 

, 2 1 579E+02 

-.221506-01 

, 26395E +0 1 

.11332E+01 

.26933E+00 

,221576+02 

-.228866-01 

,310176+01 

,139616+01 

,300776+00 

.22661E+02 

-.261916-01 

,337b7E+01 

.15771E+01 

,332526+00 

.230616+02 

-.256196-01 

,366596+01 

,182596+01 

.36960E+00 

,233656+02 

-.259996-01 

, 39688E+0 1 

,209076+01 

,396996+00 

.236266+02 

-, 272226-01 

.92876E+01 

.23720E+01 

,029716+00 

,238326+02 

-.30539E-01 

,962276+01 

,266866+01 

,062766+00 

,239936+02 

-.312796-01 

,997366+01 

,298006+01 

.09613E+00 

.28116E+02 

-.329756-01 

,539 l 3E+0 1 

,330556+01 

,5298«E+00 

,202066+02 

-.335666-01 

.57260E+01 

,369<j3E + oi 

.563896+00 

,29 2686+02 

-.396216-0) 

, 6 1 28 1 E+0 1 

,399606+01 

.59827E+00 

,293076+02 

-.356906-01 

.65979E+01 

,935986+01 

,633006+00 

.29326E+02 

-.366266-01 

, 69857E+0 1 

,973516+01 

,668086+00 

,293286+02 

-.375786-01 

,799166+01 

.51219E+01 

,703516+00 

,2931 6E+02 

-, 38999E-0 1 

.79159E+01 

, 55 1 80E + 0 1 

,739306+00 

.29291E+02 

-.39375E-01 

,890876+01 

.59295E+01 

,775996+00 

,292566+02 

-.902186-0 1 

,892096+01 

.63902E+01 

,811996+00 

,292136+02 

-.91 022E-0 1 

,995116+01 

.676976+01 

,898816+00 

.291626+02 

-.917866-01 

,100016+02 

, 7 1 97*16 + 0 1 

,886096+00 

,291096+02 

-.925086-01 

,105706+02 

.76378E+01 

,923656+00 

.29092E+02 

-.931886-01 

,11 1 58E + 02 

,808556+01 

,961696+00 

,239756+02 

-, 93823E-0 1 

,117676+02 

,853996+01 

, 10000E + 01 

.239096+02 

-.999196-01 

.12395E+02 

.900076+01 
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TABLE I.- Continued 


(c) Sphere; K = 0.823 


n 


o. 

.2879BE-01 
.S7789E-01 
,871 10E-01 
.11673E+00 
.196698+00 
,176868+00 
.207388+00 
,2 3620E+00 
.26933E+00 
.30077E+00 
.33252E+00 
.36960E+00 
,396998+00 
.92971E+00 
.96276E+00 
, 996 1 3E + 00 
.529898+00 
,563898+00 
,59827fc+00 
,633008+00 
.666088+00 
,703518+00 
,739308+00 
,775998+00 
,81 199E+00 
,898818+00 
.88609E+00 
.92365E+00 
.96169E+00 
,100008+01 


h l,1 


,800028+01 
,800038+01 
, 800 09E+0 1 
,800098+01 
,800058+01 
.60005E+01 
.80005E+01 
.80005E+01 
,800068+01 
,800068+01 
,800068+01 
,800068+01 
, 800068+01 
,800068+01 
, 8 0 0 0 & E + 0 1 
,800068+01 
,800068 + 0 1 
,800068+01 
.B0006E+01 
. 80006 E +01 
.80006E+01 
,800078+01 
,800078+01 
.80007E+01 
.80007E+01 
.80007E+01 
.80007E+01 
.80007E+01 
.80007E+01 
,800078+01 
.80007E+01 


h l,2 


,265198-02 
, 29355E-02 
, 32236E-02 
,351558-02 
,381088-02 
,910998-02 
,9911 IE-02 
.97156E-02 
,502318-02 
,533328-02 
.56961E-02 
,596168-02 
.62797E-02 
,660058-02 
.69238E-02 
,729978-02 
.757818-02 
,790928-02 
.82928E-02 
.857908-02 
.69176E-02 
,925928-02 

.96032E-02 
, 99999E-02 
, 1 0299E-0 1 
,106518-01 
,110068-01 
.113638-01 
, 1 1723E-0 1 
, 12086E-01 
» 12952E-0 1 


h 2, 2 


, 9997 1 E + 00 

.117178+01 

.135928+01 

, 15626E+01 

.178298+01 

,201 89E+0 1 

.22727E+01 

.25993E + 0 1 

.28390E+01 

,319258+01 

.39701E+01 

.38179E+01 

.916508+01 

, 95733E + 0 1 

.99630E+01 

.59195E+01 

,586858+01 

.639598+01 

.68961E+01 

.73709E+01 

.79207E+01 

.899S9E+01 

.90973E+01 

,972558+01 

.10381E+02 

. 110658+02 

,117788+02 

.125218+02 

, 13299E+02 

.19098E+02 

,199398+02 


h 3, 3 


,962698+00 
.S90858+00 
,626098+00 
.718398+00 
,8181 0E+00 
.92538E+00 
,109098+01 
,1163«8+01 
,129978+01 
.19393E+01 
,158268+01 
.173988+01 
,190618+01 
.20817E+01 
.22669E+01 
.29620E+01 
.26672E+01 
,28626E+01 
, 31091 E+ 01 
,3396*E+01 
, 35996E+0 1 
,385998+01 
,912618+01 
, 99099E + 0 1 
,970608+01 
.50199E+01 
, 53366&+0 1 
.56721E+01 
,6021 1E+01 
, 6389 1 E + 0 1 
,6761 6E + 0 1 
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TABLE I.- Continued 


(d) Sphere; K - 2.07 


n 

h l,l 

h l ,2 

h 2,2 

h 3,3 

o. 

,8001 5E+0 1 

.,192336-02 

, 99955E + 00 

.81905E+00 

.28798E-01 

,800156+01 

- , 26339E-02 

, 1 1 897E + 0 1 

.952176+00 

,577896-01 

,800 15E+01 

-.33281E-02 

, 13B72E + 01 

.110266+01 

.871106-01 

,800 15E+0 1 

-.90103E-02 

, 1 6075E + 0 1 

,126586+01 

, 1 1673E+00 

.80015E+01 

• 1 46825E»02 

, 18961E + 01 

.199206+01 

, 19669E+00 

,800156+01 

-.53966E-02 

,210396+01 

,163186+01 

,176866+00 

.80015E+01 

-.60091E-02 

.23799E+01 

,183506+01 

.20738E+00 

.80015E+01 

-,665606-02 

,26760E+01 

,205296+01 

,238206+00 

.80015E+01 

-.730316-02 

.299226+01 

.22895E+01 

,269336+00 

.80015E+01 

-,799606-02 

,332916+01 

,2531“E+01 

,300776+00 

,800156+01 

-,858536-02 

,368706+01 

.27936E+01 

, 332526+00 

,800156+01 

-.92219E-02 

, 90666E + 0 1 

,307186+01 

,36«60E+00 

,800156+01 

-.965956-02 

,996836+01 

.336576+01 

,396996+00 

,800156+01 

-, 10985E-01 

.98928E+01 

.367696+01 

,9297 16+00 

,800196+01 

-.111136-01 

,53905E+01 

,900916+01 

.96276E+00 

.60010E+01 

-.117396-01 

,581206+01 

,939936+01 

.99613E+00 

.80019E+01 

• , 1 2363E-0 1 

,630806+01 

.97129E+01 

.52989E+00 

,800196+01 

• , 1 2985E-0 1 

,682896+01 

,5093»E+01 

,563896+00 

.80019E+01 

- , 1 3605E»0 1 

,737596+01 

,599916+01 

,598276+00 

,800196+01 

-.19223E-01 

,799826+01 

.591386+01 

.63300E+00 

,800 1 9E+0 1 

- , 1 9B39E-0 1 

,859786+01 

,635396+01 

,668086+00 

,8001 9E+0 1 

-.15959E-01 

,917996+01 

,681336+01 

,703516+00 

.80019E+O1 

-, 160676-01 

,983016+01 

.72991E+01 

.73930E+00 

.80019E+01 

-, 16679E-01 

,1051 9E+02 

,779636+01 

,775996+00 

,800196+01 

-,172896-01 

,112286+02 

,832056+01 

,81 199E + 0Q 

,800196+01 

-.17898E-01 

, 1 1972E+02 

.88672E+01 

,898816+00 

,800 19E+01 

-.185056-01 

, 12796E+02 

,993716+01 

, 88609E+00 

,8001 9E+0 1 

-.191116-01 

,135536+02 

.10031E+02 

, 92365E + 00 

.80019E+01 

-.19715E-01 

,193926+02 

.10698E+02 

,961696+00 

.80019E+01 

-.203176-01 

,152696+02 

.U291E + 02 

, 10000E+01 

,800196+01 

-.209186-01 

.161706+02 

.U959E+02 
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TABLE I.- Continued 


(e) Blunt body; K - 0.791 


n 

h l,l 

h l ,2 

h 2, 2 

^3,3 

o. 

.987226+01 

,260706-03 

,880076+00 

.440816+01 

.287486-01 

,988706+01 

,235736-03 

,916666+00 

,474256+01 

.57784E-01 

.98962E+01 

,186836-03 

.95726E+00 

.S09496+01 

.871 10E-01 

,990096+01 

,159856-03 

,100186+01 

.546586+01 

, 1 16736+00 

.99021E+01 

,141656-03 

.10501E+01 

.585576+01 

, 146646+00 

,960066+01 

,131 35E-03 

. 1 1022E+01 

,626496+01 

,176866+00 

,98967E+«l 

.127346-03 

,115816+01 

, 6694 l E + 0 1 

.207386+00 

, 9891 OE + 0 1 

.129036-03 

, 121 77E+01 

,714376+01 

.236206+00 

.98837E+01 

.135786-03 

, 12812E+01 

,761436+01 

.269336+00 

.98751E+01 

, 1 6704E-03 

, 13485E+01 

,810656+01 

,300776+00 

.98655E+01 

,162286-03 

, 1«19#E+01 

,862076+01 

,332526+00 

,985496+01 

,180946-03 

,149516+01 

.915786+01 

,364606+00 

.98436E+01 

,202«5E-0 3 

, 15745E+01 

,97 179E+0 1 

,396996*00 

,903156+01 

, 22620E-03 

.16S80E+01 

.10302E+02 

,429716+00 

,981896+01 

,254866-03 

, 17458E+01 

,1091 1E+02 

,<<62766 + 00 

,980586+01 

, 28456E-03 

, 1 8379E+0 1 

.11545E+02 

, <496135 + 00 

,979236+01 

,484426-03 

,193456+01 

, 12204E+02 

,529846+00 

.97784E+01 

.52308E-03 

.20359E+01 

.128916+02 

.56389E+00 

.97641E+01 

,556146-03 

,214206+01 

.13604E+02 

,598276+00 

.97496E+01 

,581606-03 

•22528E+0 1 

, 143«8E+02 

,633006+00 

.97348E+01 

,537976-03 

,236876+01 

,151186+02 

,668086+00 

,971986+01 

-.984726-03 

, 24895E+0 1 

.15915E+02 

,703516+00 

.97046E+0! 

-.61 374E-03 

,261576+01 

.18745E+02 

.73930E+00 

.96H92E+01 

-.720916-03 

.27471E+01 

.17605E+02 

,775446+00 

.96755E+01 

-.795066-03 

, 28840E + 0 1 

,180976+02 

,811946+00 

.964516+01 

-.57918E-03 

.30263E+01 

,194226+02 

,848816+00 

,962836+01 

-.36024E-03 

.31743E+01 

.20379E+02 

,886046+00 

.96271E+01 

-,384996-03 

,332746+01 

,213706+02 

,923656+00 

,960976+01 

-.379876-03 

,348636+01 

.22397E+02 

.961646+00 

,959416+01 

-.35825E-03 

,365116+01 

,234606+02 

,100006+01 

.95785E+01 

-.32303E-03 

,3821 9E+0 1 

,245596+02 
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TABLE I.- Concluded 


(f) Blunt body; K - 4.28 


n 


0 « 

.28748E-01 

.S7784E-01 

.871 10E-01 

,ilb73E+00 

. 14bb4E+00 

, 17686E+00 

.20738E+00 

.23820E+00 

,2b933E+00 

,30077Ef0u 

, 33252E+00 

>,3b4b0t«-00 

.39699E+00 

.42971E+0Q 

,«6276E+00 

,49b 1 3E + 0 0 

,52984E+00 

,Sb389£+00 

.59827E+O0 

.83300E+00 

,bb808E+00 

.70351E+00 

.73930E+00 

,77544Et00 

.81194E+00 

, 8488 1E + 00 

,B8b04Ef00 

,923b5E+Oo 

,9t»lb4t + 00 

, 10000E + 01 


h l,l 


( 1 1 203E+02 
, 1 1 199E+02 
,11 1 93E + 02 
.11195E4-02 
, 1 1 194E+02 
,11 1 8BE+02 
, 1 1 187E + 02 
, 1 1 1 8bE+02 
,111 8bE+02 
.1118SE+02 
, 1 1 185E+02 
,11 185E+02 
,111 85E+02 
,11 185E+02 
,111 8bE+02 
,11 181E+02 
, 1 1221E + 02 
, 1 1 22 1 E + 02 
, 1 1 1B0E+O2 
,1 1 184E + 02 
, 1 1 183E + 02 
, 1 1 1 82E + 02 
,11 181E+02 
,11180Et02 
,11 180E+02 
,11 179E+02 
, 1 1 179E + 02 
,111 78E+02 
.11 178E + 02 
,11 178E+02 
,11 177E + 02 


h l ,2 


,2b413E-02 

,285b0E-02 

.33395E-02 

,37b38E-02 

,b7931E-02 

,32l5bE-02 

.49259E-02 

,4738bE-02 

,4b72«E-02 

,4b309E-02 

,4b2b8E-02 

,4b710E-02 

.47717E-02 

.49350E-02 

.52830E-02 

.78549E-02 

.78449E-02 

.80243E-02 

.85494E-02 

,b2bbbE-02 

,b8510E-02 

,b8398E-02 

.89349E-02 

.70817E-02 

.77385E-02 

.78999E-02 

.80783E-02 

.82904E-02 

.8S405E-02 

.B8303E-02 

.91693E-02 


h 2,2 


, 10913E+01 
, 12241E+01 
, 1 3649E + 0 1 
.15138E+01 
• 1 67 1 1E + 0 1 
.18362E+01 
,201 0 0E + 0 1 
.21917E+01 
,23821Et01 
.25812E+01 
.27892E+01 
,300b4E+01 
, 32330E+0 1 
,34b92Et01 
, 37 1 53E+0 1 
.39721E+01 
. 42393E + 0 1 
.95184E+01 
.48100E+01 
.S1137E+01 
, 54287E + 0 1 
, 57556E+0 1 
,b0947E+01 
, b44b3E + 0 1 
, b8 1 07E + 0 1 
.71878E+01 
.75784E+01 
, 79828E + 0 1 
,840 1 bE+0 1 
, 8835 1 E+0 l 
.92839E+01 


h 3 f 3 


,14889Et01 
,lbl67Ef0i 
,1752/E+Ol 
,189b8£+01 
.20491E+01 
.22098E+01 
.23793E+01 
.25579E+01 
,27458Ef01 
.29433E+01 
.31507E+01 
, 33883E+0 l 
, 35985E+0 1 
.3835&E+01 
,4085bE+01 
.43473E+01 
,4620bE+0i 
,4906®E+0 1 
.52046E+01 
,5515bE+01 
, S8397E + 0 1 
,bl779E+01 
, 85290E + 0 1 
.68950E+01 
,7275«E+01 
,7b7l7E+01 
, 80833E + 0 1 
.85109E+01 
, 89550E + 0 1 
,94ifelEt01 
, 9894&E + 0 1 







TABLE II.- COORDINATE LINE SLOPES 


(a) Ellipse 


n 

K - l 

.34 

K = 3 

.72 

(dy/dx) n 

/ (dy/dx) ^ 

(dy/dx) n 

-V (dy/dx) j- 

.287486-01 

.20326E+0O 

,19‘924EtOO 

,91493EtOO 

.89472E+00 

,577846-01 

,232l6Et00 

,22923Et00 

,91737Et00 

,901 05E + 0 0 

.87U0E-01 

.256836+00 

,25473Etoo 

,93828Et00 

,923486+00 

, 116736+00 

.27847E+00 

.276846+00 

,96974EtOO 

.95540E+00 

, 146646+00 

,29768Et00 

,29632Et0O 

, 1 0069EtO 1 

,993176+00 

.176866+00 

.31519E+00 

,31387Et00 

, 10481Et01 

, 1 0346E + 0 1 

,207386+00 

• 331 29E+00 

, 3 3 0 1 2 E 1 0 0 

.109UE + 01 

, 1 0782E + 0 1 

,238206+00 

.345B6E+00 

, 345 1 OEtOO 

,1 1360E + 01 

, 1 1231 E + 0 1 

.269336+00 

,35988EtOO 

,3591 1 E tOO 

,118196+01 

,1 1690E+01 

.30077E+00 

,37314E+0O 

,37244£t00 

, 1 2284E + 0 1 

. 121 59 E + 01 

.33252E+C0 

,38577Et00 

, 3851 3Et00 

.12763E+01 

, 12631E + 0 1 

,364606+00 

,397«8Et00 

,39728Et00 

, 1 3236E+0 1 

,131006+01 

,396996+00 

,40955£t00 

,40899EtOO 

.137116+01 

, 13570E + 0 1 

,4297 16 + 00 

, 42085 E tOO 

,420326+00 

.14182E+U1 

.14041E+01 

.46276E+00 

,43l84EtOO 

,43l33Et00 

, 14662E+01 

,145146+01 

,496 1 3E + 0 0 

,44257Et00 

,44208Et0O 

,151 45E + 0 1 

, 1 499 1 E + 0 1 

.52984E+00 

,45308Et0D 

,45259EtOO 

, 1 5632E + 0 1 

, 15470E + 0 1 

,563896+00 

,U63«OE+00 

,46292Et00 

,161236+01 

, 15953E + 0 1 

.598276+00 

,47355Et00 

, 47308Et00 

, 16620E + 01 

, 1644 1E + 0 1 

,633006+00 

,48358EtOO 

,4831 lEtOO 

,171216+01 

.16933E+01 

,668086+00 

,49348Et00 

,49302Et00 

.17629E+01 

,174306+01 

,703516. tOO 

,50329Et00 

,5O283Et00 

, 18143E+01 

,179336+01 

,73930Et00 

,513026+00 

,51255Et00 

, 1 8664E + 0 1 

,184436+01 

,77544Et00 

,52269EtOO 

,52221 EtOO 

.191936+01 

, 1 8960E+0 l 

,811946+00 

.53230E+O0 

,5318lEt00 

, 1 9730E + 0 1 

,194846+01 

,848816+00 

,541 BbEtOO 

,54136EtO0 

, 20276E+0 1 

,200176+01 

,886046t00 

,55139Et00 

.55087E+00 

,208326+01 

.2O558E+01 

,92365Et00 

,56089EtO0 

,560356+00 

,21 398E+0 1 

,211086+01 

,96164Et00 

,5709SEt00 

,56995Et00 

i 

.219746+01 

.21668E+01 
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TABLE II.- Continued 


(b) Sphere 


n 

K = l 

D.823 

K = 2.07 

(dy/dx) n 

/ (dy/dx) ^ 

(dy/dx) n 

-1/( dy/dx) £■ 

.28748E-01 

.57784E-01 

.87110E-01 

,116736+00 

* 146646+00 

, 1 76866+00 

, 20738E+00 

,23820E+00 

.26933E+00 

,300776+00 

.33252E+00 

.36460E+00 

.396996+00 

#4297 1E + 00 

.462766+00 

.496138+00 

.529848+00 

, 56389E + 00 

.59827E+00 

.63300E+00 

,668086+00 

. 70351E + 00 

.73930E+00 

.77544E+00 

,811946+00 

.84881E+00 

.88604E+00 

. 9236SE+00 

.961646+00 

. 1 0821E+01 
. 1 0823E+O1 
. 1 0H25E + 0 1 
.108266+01 
. 10828E+01 
.108306+01 
. 10831E+01 
. 1 0832E + 0 1 
. 1 0834E + 0 1 
■ 1 0835E + 0 1 
. 1 0836E+0 1 
.10837E+01 
. 10838E+01 
. 1 0839E + 0 1 
. 10840E+01 
.10841E+01 
« 1 0842E + 0 1 
.10843E+01 
. 1 0844E + 0 1 
, ) 0B45E+0 1 
, 1 0846E + 0 1 
. 1 0847E + 0 1 
. 108O8E + 01 
, 10849E + O1 
, 1 0849E + 0 1 
. 1 0850E + 0 1 
.108516+01 
. 1 0852E + 0 1 
.108526+01 

.107996+01 
.10801E+01 
. 1 0802E + 0 1 
. 1 0804E + 0 1 
, 10805E + 01 
, 10807E + 01 
,108086+01 
,108096+01 
,108116*01 
,108126+01 
.108136+01 
,108146+01 
,108156+01 
,108166+01 
,108176+01 
,108186+01 
, 10819E + 01 
,108206+01 
, 10821E + 0 1 
,108226+01 
.108236+01 
,108236+01 
,108246+01 
,108256+01 
, 10826E + 01 
,108276+01 
, 1 0827E + 01 
,108286+01 
,108296+01 

, 1 8585E + 0 1 
, 18583E+01 
,185826+01 
. 18580E+01 
, 1 8578E + 0 1 
, 1 8577E+0 1 
,185756+01 
.185746+01 
, 1 8573E + 0 1 
, 1 857 1 E + 0 1 
.185706+01 
, 18568E+01 
, 18567E+01 
,185666+01 
.185656+01 
. 16563E + 01 
.185626+01 
. 18561E + 01 
. 185606 + 01 
. 1 8558E + 0 1 
,185576+01 
.185566+01 
.185556+01 
• 1 8554E + 0 1 
, 1B553E+01 
. 1 8552E+0 1 
. 18551E + 01 
.185506+01 
, 18548E + 01 

.184226+01 
.184196+01 
.184J6E+01 
,184136+01 
,184106+01 
.184076+01 
.18405E+01 
,184026+01 
,184006+01 
. 18398E + 01 
,183966+01 
.183946+01 
.183926+01 
.183906+01 
.183886+01 
,183866+01 
, 183856 + 01 
.183836+01 
,183816+01 
,183806+01 
, 18378E + 01 
.183776+01 
. 18375E + 0 1 
.18374E+01 
, 1 8372E+0 1 
.183716+01 
.18370E+01 
, 1 8368E + 0 1 
.18367E+01 
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TABLE II.- Concluded 


(c) Blunt body 



K = c 

1.791 

C - 4. 

28 

n 

(dy/dx) n 

-1 / (dy/dx) ^ 

(dy/dx) n 

-1 / (dy/dx) ^ 

,287486-01 

.572416+00 

.57270E+00 

, 15944E+01 

,156956+01 

.577846-01 

,559536+00 

,559776+00 

.157596+01 

, 15521 E+0 1 

,871106-01 

.54757E+00 

,547766+00 

.155906+01 

.153656+01 

, 1 1673E+00 

.53637E+00 

.536526+00 

,154596+01 

. 15234E+0 1 

, 146646+00 

.52579E+00 

,525916+00 

.15290E+01 

.15107E+O1 

, 1768b6+00 

.51574E+00 

,515836+00 

.151586+01 

,149786+01 

,207386+00 

.506156+00 

.50622E+00 

. 15036E+U1 

.148616+01 

,238206+00 

,4969bE+00 

.49701E+00 

,149226+01 

.147526+01 

,2 69336+00 

.488126+00 

,488146+00 

.148156+01 

,146496+01 

,300776+00 

,479586+00 

,479596+00 

.14715E+01 

, 14552E + 0 1 

, 332526+00 

.471316+00 

,471316+00 

.146216+01 

, 14460E + 0 1 

,364806+00 

,463296+00 

,463286+00 

.14532E+01 

,143736+01 

,396996+00 

,455506+00 

,455476+00 

, 14448E + 01 

.14291E+01 

,429716+00 

,447916+00 

,447876+00 

, 1 4368E + 0 1 

,1421 3E+01 

.46276E+00 

,440506+00 

,440456+00 

.143006+01 

.141416+01 

,496136+00 

.43327E+00 

,43321 E + 00 

.142286+01 

,140746+01 

.52984E+00 

,426206+00 

,426136+00 

. 1 4165E + 01 

.140076+01 

,563896+00 

.41928E+00 

,419216+00 

,140986+01 

.13941E+01 

,598276+00 

.412516+00 

,412426+00 

.140236+01 

. 138766 + 01 

,633006+00 

.405B7E+00 

.40576E+00 

. 13960E + 01 

,138146+01 

,668086+00 

.39970E+00 

.39933E+00 

,139006+01 

,137556+01 

, 70351E+00 

. 39332E + 00 

.39301E+00 

.1384)6+01 

.136986+01 

.73930E+00 

.3B704E+00 

.38673E+00 

. 13784E + 01 

.136436+01 

,775446+00 

,380886+00 

,380556+00 

.137296+01 

,135906+01 

.81194E+00 

.37481E+00 

.374516+00 

.136766+01 

. 1 3539E + 0 1 

,8488lE+00 

.36884E+00 

, 36858E+00 

, 1 3624E + 0 1 

,134896+01 

,886046+00 

,362966+00 

.36269E+00 

, 1 3574E + 0 1 

.134406+01 

,923656+00 

.357176+00 

.35689E+00 

.135256+01 

, 1 3392E + 0 1 

,961646+00 

.351476+00 

,351196+00 

, 1 3478E + 0 l 

,133466+01 
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TABLE III.- METRIC COEFFICIENT AS FUNCTION OF MESH SIZE 

FOR 1.5:1 ELLIPSE 



^1,1 for coordinate mesh size of - 


101 x 41 

101 x 31 

101 X 21 

51 x 41 

51 x 31 

51 x 21 

0 

13.000 

13.000 

13.000 

13.000 

13.000 

13.000 

.159 

12.996 

12.996 

12.996 

12.996 

12.996 

12.996 

.317 

12.567 

12.568 

12.570 

12.581 

12.581 

12.583 

.456 

12.349 

12.349 

12.350 

12.190 

12.190 

12.192 

.635 

11 .842 

11 .843 

11 .846 

11 .750 

11 .750 

11 .753 

.793 

11 .344 

11 .345 

11 .349 

11.247 

11.248 

11 .251 

.872 

10.796 

10.797 

10.831 

10.668 

10.669 

10.673 

.912 

10.302 

10.303 

10.324 

10.137 

10.138 

10.142 

1 .031 

9.857 

9.858 

9.862 

9.686 

9.686 

9.690 

1 .189 

9.555 

9.556 

9 . 597 

9.366 

9.366 

9.369 

1 .348 

9.414 

9.414 

9.432 

9.239 

9.239 

9.241 

1 .507 

9.453 

9.453 

9.454 

9.315 

9.316 

9.317 

1 .666 

9.583 

9.583 

i i 

9.584 

9.524 

9.524 

9.525 


2 
















Figure 1 Body surface geometry. 


21 




Figure 2.- Geometry of body and outer boundary. 
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Figure 5.- Orthogonal coordinate system abput 1.5:1 ellipsoid with 

equal £ spacing. 


Figure 6.- Orthogonal coordinate system about 2:1 ellipsoid with 

equal £ spacing . 


Figure 7.- Orthogonal coordinate system about 3:1 ellipsoid with 

equal % spacing. 



Figure 8.- Orthogonal coordinate system about 3:1 ellipsoid with slightly 

unequal £ spacing. 
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Figure 9.- Orthogonal coordinate system about 3:1 ellipsoid with 

unequal £ spacing. 
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Figure 10.- Orthogonal coordinate system about 3:1 ellipsoid with 

unequal K and n spacing. 
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Figure 11.- Orthogonal coordinate system about planetary entry body. 


Figure 12.- Orthogonal coordinate system about planetary entry body with 

unequal n spacing . 



Figure 13.- Orthogonal coordinate system about arbitrary body with 

reverse curvature. 
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Figure 14 


Orthogonal coordinate system about arbitrary body with reverse 
curvature and unequal spacing in n direction. 
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Figure 15.- Orthogonal coordinate system about 3:1 ellipsoid with displaced 
outer boundary and unequal £ spacing. 
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Figure 16.- Orthogonal coordinate system about planetary entry body with 
displaced outer boundary and unequal r| spacing. 
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Figure 18.— Orthogonal coordinate system on forebody of reverse curvature 

conical body with unequal n spacing. 
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Figure 19.- Coordinate systems schematic. 
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